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Abstract 



(D _ 

\J-^ \ Surfaces of experimental masses of even-even and odd-odd nuclei exhibit a 

t~^ ' sharp slope discontinuity at N=Z. This cusp (Wigner energy), reflecting an 

^^ ' additional binding in nuclei with neutrons and protons occupying the same 

shell model orbitals, is usually attributed to neutron-proton pairing correla- 
^ , tions. A method is developed to extract the Wigner term from experimental 

yf^ \ data. Both empirical arguments and shell-model calculations suggest that the 

C^ ' Wigner term can be traced back to the isospin T=0 part of nuclear interac- 

tion. Our calculations reveal the rather complex mechanism responsible for 



t^^ ■ the nuclear binding around the N=Z line. In particular, we find that the 

- Wigner term cannot be solely explained in terms of correlations between the 



neutron-proton J=l, T=0 (deuteron-like) pairs. 
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One of the most interesting avenues in nuclear structure research is the study of nuclear 
properties near the particle drip lines. Theoretically, due to their unusual shell structure 
and weak binding, nuclei with highest and lowest iV/Z ratios represent a unique challenge 
for the many-body problem. 

This Letter concerns the proton-rich border of nuclear binding [N/ Zk,V). The proton drip 
line is relatively well known experimentally, and the discovery of ^°°Sn indicates the potential 
that radioactive nuclear beams will access all bound N=Z nuclei and their bound and 
unbound neighbors. This opens up the possibility of nuclear structure studies of medium- 
mass and heavy systems with N^Z. 

There are many aspects of nuclear structure which make physics near the N=Z line 
very interesting: proton and diproton emission, isospin mixing, superallowed beta decay, 
and nuclear constraints for the astrophysical rp-process - all these are timely experimental 
and theoretical themes. A unique aspect of nuclei with N=Z is that neutrons and protons 
occupy the same shell-model orbitals. Consequently, the large spatial overlaps between 
neutron and proton single-particle wave functions are expected to enhance neutron-proton 
{np) correlations, especially the np pairing. 

At present, it is not clear what the specific experimental fingerprints of the np pairing 
are, whether the np correlations are strong enough to form a static condensate, and what 
their main building blocks are. Most of our knowledge about nuclear pairing comes from 
nuclei with a sizable neutron excess where the isospin T=l neutron-neutron {nn) and proton- 
proton (pp) pairing dominate. Now, for the first time, there is an experimental opportunity 
to explore nuclear systems in the vicinity of the N=Z line which have many valence np 
pairs; that is, to probe the interplay between the like-particle and neutron-proton (T=0,1, 
Tz=0) pairing channels. 

This novel situation calls for the generalization of established theoretical models of nu- 
clear pairing. In spite of several early attempts to extend the independent quasi-particle 
formalism to incorporate the effect of np correlations in light nuclei |I|-^ (see Ref. |^ for 
an early review), no symmetry-unrestricted calculations for np pairing, based on the quasi- 
particle theory, have been carried out. 

So far, the strongest evidence for np pairing comes from the masses of N=Z nuclei. An 
additional binding (the so-called Wigner energy) found in these nuclei manifests itself as a 
spike in the isobaric mass parabola as a function of Tz=^{N — Z) (see the review and 
Refs. quoted therein). Recently, it has been demonstrated that the T=0 np-pair correlations, 
treated within the framework of the generalized quasi-particle theory, may produce such a 
cusp at N=Z [0]. Gross estimates of the magnitude of the Wigner energy come from a large- 
scale fit to experimental binding energies with the macroscopic-microscopic approach [^,|^ 
and from the analysis of experimental masses [|10|. In this study, we present a technique to 



extract the Wigner energy directly from the experimental data, and give empirical arguments 
that this energy originates primarily from the T=0 part of the effective interaction. To obtain 
deeper insight into the structure of the Wigner term, we apply the nuclear shell model to 
nuclei from the sd and fp shells. 

In the semi-empirical mass formulae (see e.g. Ref. 0), an additional binding due to the 
np-pair correlations is usually parameterized as 

B^^f = -enp{A)7r^p + Ew, (1) 



where 7r„p = i(l — 7r„)(l — iTp) and 7r„ = (—1)'^ and vr^ = (—1)^ being the nucleon- 
number parities. The first contribution to the np pairing energy in Eq. (JlD, e„p, represents 
an additional binding due to the residual interaction between the two odd nucleons in an 
odd-odd nucleus. 

The second contribution to B^^", E^, dubbed the Wigner energy, is believed to represent 
the energy of collective np-pairing correlations. It can be decomposed into two parts: 

Ew = W{A)\N -Z\+ d{A)nr,p6Mz, (2) 



The I A^ — Z I -dependence in Eq. (^ was first introduced by Wigner |TT| in his analysis of the 



SU(4) spin-isospin symmetry of nuclear forces. In the supermultiplet approximation, there 
appears a term in the nuclear mass formula which is proportional to Tgs(Tgs + 4), where 
Tgs denotes the isospin of the ground state. Empirically, Tgs=\Tz\ for most nuclei except 



for heavy odd-odd N=Z systems |T^,[T3|. Although the experimental data indicate that the 
SU(4) symmetry is severely broken, and the masses behave according to the Tgs(Tgs -|- 1) 
dependence |T^,T^, the expression of Eq. (j^) for the Wigner energy is still very useful. In 



particular, it accounts for a non-analytic behavior of nuclear masses when an isobaric chain 
crosses the N=Z line. An additional contribution to the Wigner term, the d-term in Eq. (^, 
represents a correction for N=Z odd-odd nuclei. Theoretical justification of Eq. (0) has been 
given in terms of basic properties of effective shell-model interactions P,|l3 , and also by using 



simple arguments based on the number of valence np-pairs [TyjT^ . The estimates based on 
the number of np pairs in identical spatial orbits suggest that the ratio d/W is constant and 
equal to one [|15|. A different estimate has been given in Ref. [|1^: d/W=0.56±0.27. 



In order to extract empirical information on the magnitude of pairing correlations, one 
usually applies mass relations (indicators, filters) defined as specific combinations of nuclear 
binding energies of neighboring nuclei [jlO[- The use of these filters is based on the following 
assumptions: («) the nuclear binding energy may be decomposed into a part B which varies 
smoothly with A^ and Z (at least away from the closed shells) and a fluctuating (and possibly 
non- analytic) term; and (m) the pairing energy may be parameterized as a smooth function 
of A^ and Z to reproduce global (large-scale) behavior. The important feature of mass 
indicators, based on finite-differences of the nuclear binding energy with respect to A^ and 
Z, is that they remove the smooth energy background B up to third [|13)0 (or higher 



jlTl ) orders, and isolate specific parts of the pairing energy. However, the interpretation 
of results depends on the particular indicator used and the assumed parameterization of 
pairing energy. 

The magnitude of the Wigner term is non-negligible only in the closest vicinity of the 
N=Z line. However, since pairing indicators are given by finite differences, perturbations 
at N=Z might disturb the physical interpretation of extracted quantities. For example, the 
well-known four-point odd-even mass difference formulae for the nn and pp pairing gaps are 
quite strongly perturbed near the N=Z line |]18|. Consequently, to isolate interesting effects. 



other indicators, involving more nuclei, must be used. 

An indicator particularly useful for investigating the np correlation energy is the double- 
difference formula of Ref. |T^ : 

SVnp{N, Z) = - {B{N, Z) - B{N -2,Z)- B{N, Z - 2) + B{N -2,Z - 2)} 



dNdZ' 



(3) 



This indicator belongs to a different class than those analyzed in Ref. |]T0|- First, it in- 
volves binding energies of nuclei having the same nucleon-number parity and, therefore, it 
is sensitive to the Wigner energy |20|. Second, indicator (|^) does not remove the second- 



order contribution to the smooth energy proportional to NZ . As a consequence, for nuclei 
with Tz > 1, SVnp mainly probes the symmetry energy term, ~ (A^ — Z)'^/A, and for nu- 
clei with T^ < 1 it contains both Wigner energy and symmetry energy contributions. The 
nuclear symmetry energy is well described by a self-consistent mean-field theory based on 
the realistic effective particle-hole interaction. Indeed, the application of indicator (§) to 
the theoretical self-consistent mass tables using either a spherical Hartree-Fock-Bogoliubov 
model (HFB) ||2l| or a deformed extended Thomas-Fermi plus Strutinsky-integral model 

for nuclei 



2| gives results very consistent with the experimental values of 6Vr 
1 



np 



T8[] . This simple test indicates that the symmetry energy is rather insensi- 



(ETFSI) 

with Tz > 

tive to shell effects. Its strength smoothly varies with A and is practically T^-independent. 

Since, in this study, we are mainly interested in the Wigner energy, another indicator has 

to be employed that allows for filtering out the symmetry energy contribution. This can be 

achieved by considering combinations of SVnp values, as discussed in the following. 



In this work, the Wigner energy coefficient W in an even-even nucleus Z=N= 
extracted by means of the indicator: 



W{A) = 6Vr 
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The (i-term in an odd-odd nucleus, Z 
cator: 
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[Eq. (0)] can be extracted using another indi 
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(5) 



Although the recipe for these third-order mass difference indicators is not unique, the results 
appear to be very weakly dependent on the particular prescription used |]18 |. 



The experimental values of both contributions to the Wigner energy are shown in Fig. |l|. 
The values of W decrease rather smoothly with A, showing characteristic oscillations that 
are related to the shell effects. Indeed, the filling of l(i5/22si/2 shells results in a decrease 
of W with the local minimum near y4~34. Subsequent filling of the 1^3/2 shell tends to 
increase W. The next minimum in W appears at the middle of the I/7/2 shell. Table | 
shows the results of a least-squares fit to W assuming Vr(yl) = aw /A"'- The actual minimum 
corresponds to q;~0.95; it is very close to the standard value of a=l. The extracted strength 
avi/~47MeV is considerably larger than the one used in the semi-empirical mass formula of 
Ref. (avK=30MeV), but is consistent with the value of Ref. |T0| (a^=43±5MeV). 

As seen in Fig. |l|, for the sd nuclei, the behavior of the (i-term follows rather closely 
that of W, but this nice agreement is lost for the heavier systems. However, it is known 
experimentally (cf. discussion in Ref. and references quoted therein) that as a rule the 
nuclear ground states have isospin Tgs=|A^ — Z\/2, except for heavy N=Z odd-odd nuclei. 
Indeed, ^"^Cl and nuclei heavier than ^°Ca (except for ^^Cu) have ground states of isospin 
T=l. Replacing, in our analysis, the ground-state binding energies of these odd-odd N=Z, 



T=l nuclei by the binding energies of their lowest T=0 states leads to modified (i- values 
marked by triangles in Fig. |l]. The resulting values follow very closely the values of W; 
i.e. dT=o{A)/W{A) !=a 1, in a nice agreement with the simple arguments of Ref. [|T5|. This 
observation provides an important experimental argument that both terms constituting E\y 
have the same microscopic origin. In particular, assuming that the structure of the lowest 
T=0 states in odd-odd N=Z nuclei is strongly infiuenced by the effective interaction between 
T=0 pairs, the close similarity between dT=o and iy(y4) indicates that the T=0 component 
of the nuclear interaction is responsible for the Wigner term. 

The role of the T=0 part of the interaction on the presence of binding-energy irregularities 
in 6Vnp near the N=Z has been recognized in Ref. |^. They nicely illuminated this point 
by demonstrating that in the shell-model calculations for sd nuclei using the surface delta 
interaction (SDI), it is indeed the T=0 part of SDI that give rise to the cusp in SVnp. 

In order to investigate the role played by the T=0 component for the Wigner energy, and 
to study the relative importance of various np pairs, we performed shell-model calculations 



in both the lf2p and ld2s shell-model spaces using the shell-model code ANTOINE |24]. In 



the /p-shell we use the KB3 effective interaction of Ref. [^ , while for the sd-shell nuclei we 



use the USD effective interaction of Ref. ||2^. The results of the calculations are displayed 
in Figs. and ^ 

The insert in Fig. Q shows the calculated binding energies (solid circles) for the A=48 
isobaric chain of even-even fp nuclei relative to the N=Z nucleus ^^Cr. The Wigner term 
gives rise to an additional binding for the N=Z nuclei. In Fig. 0, this is represented as a 
deviation, ew, from a smooth parabola fitted to the N^Z isobars. Note, however, that the 
quantity, ew, is not a realistic estimate of the Wigner energy Ew, but a rather schematic 
measure of the cuspQ 

To visualize the infiuence of the T=0 part of the effective nuclear interaction on the 
Wigner term, we have performed a set of shell-model calculations while switching off se- 
quentially the J = 1,2, ..., Jmax, T=0 two-body matrix elements {jiJ2JT\H\j[J2JT) of the 
shell-model Hamiltonian H for different values of Jmax- Figure ^ shows a ratio ew/^w'^^^ 
where e]^^^ denotes the result of full shell-model calculations versus Jmax- The calcula- 
tions were performed for two representative examples, namely, the fp-sheW nucleus "^^Cr and 
the sd-shell nucleus ^^Mg. The largest contribution to the Wigner energy comes from the 
part of the T=0 interaction between deuteron-like (J=l) and 'stretched' [J=5 (sd) and 7 
(pf)] pairs. The importance of these matrix elements is well known; it is precisely for J=l 
and stretched pair-states that experimentally determined effective np T=0 interactions are 



strongest [p^ -p9[|. Note also that the deuteron-like correlations contribute more strongly 
to Ew in sd-nuclei than in /p-nuclei, and that matrix elements corresponding to interme- 
diate values of J give non-negligible contributions. This reveals the complex structure of 
the Wigner energy, and suggests that models which ignore high-J components of the np 
interaction (e.g., by considering only J=0, T=l and J=l, T=0 np pairs ||30|) are not too 



^In reality, there is some contribution from the collective np-pairing effects also in nuclei with 
|A^ — Z\>4:, and, due to the isospin invariance of the interaction, the isotopic chains of even-even 
nuclei considered in Fig. |2| contain only three independent masses. Hence, the simple parabolic fit 
does not represent the realistic symmetry-energy dependence. 



useful for discussing the actual np pair correlations. 

Since in the self-consistent models based on the realistic effective interactions the explicit 
np coupling is ignored, they cannot reproduce the Wigner term (see, however, Ref. [0; 
for the recent extension of the HFB method to np pairing, see Ref. |3ll]). The successful 
ETFSI mass formula, based on Skyrme-like particle-hole interaction and monopole-pairing 
interaction between like nucleons, accounts well for most effects attributed to the particle- 
hole channel, e.g., shell effects and deformations. However, all those are filtered out by 
indicator (^. This is illustrated in the insert to Fig. ^ which shows that iy~0 in ETFSI. 

The Qhuj shell model, on the other hand, reproduces well experimental values of W both 
in the sd- as well as in the /p-shell (circles in Fig. ^. The results of calculations with the 
( J=l, T=0) two-body matrix elements removed, marked by triangles in Fig. §, demonstrate 
that the impact of the deuteron-like correlations on the strength of the Wigner term is 
stronger in sc?- nuclei than in /p-nuclei (cf. also Fig. ^). We also note the strong variations 
in the region of y4~32 where W becomes locally negative. These variations are associated 
with the subsequent filling of the upper part of the 1^5/2, the 2si/2, and of the bottom part of 
the l(i3/2 shells. In the extreme shell-model picture, both 1^5/2 and 2si/2 shells are occupied 
for Z=N=1Q, and the ^3/2 shell is well separated. Consequently, the important contribution 
to W comes from T=0 matrix elements with jn^jp)- 

Removing all T=0 matrix elements with J=l,3,5,7 {J^^=7 variant), i.e. those which 
are predominantly due to the coupling between neutrons and protons in identical orbits, 
washes out W almost completely. An exception is the upper part of the sii-shell (1^3/2 
subshell) where the value of W is reduced by about 50-60%. Intuitively, below ^°Ca one can 
expect a quenching of phase-space available for scattering the np-pairs in identical orbits. 
Simultaneously, one would expect a slight increase of the role of jn^jp rip-pairs. Indeed, 
removing the remaining T=0 matrix elements with even values of J (Jmax=7 variant) has 
the largest impact on the heaviest srf-shell nuclei. 

We briefly comment on the first term, e„p(y4), in Eq. (|I|). The magnitude of the residual 
np correlation between the last nucleons in an odd-odd system may be extracted from 
experimental masses using the pairing indicators ^np(pn) proposed in Ref. 0] . In the semi- 
empirical mass formulae, the most common parameterization of enp{A) assumes a volume- 
like dependence of enp{A) ^ e/A with eRi20-40MeV [pjp^. Our analysis performed for 
A>2Q shows that the enp{A) is almost mass-independent for nuclei around the N=Z line 
(with \N — Z\ < 3), and is of the order of e np{ A) ^500 keV . Therefore, particularly for light 
nuclei, the magnitude of Enp is 2-3 times smaller than suggested in Refs. PJIDj]. For heavier 
nuclei with larger neutron excess, Enp is systematically weaker and approaches the value of 
-300 keV. 

In summary, we performed a systematic study of rap-pair correlations in Nr^Z nuclei. 
Our main emphasis has been the analysis of the Wigner term in terms of np pairs. The 
experimental data are consistent with a two-term structure given by Eq. (^). We developed 
a reliable method which allows for the extraction of both W and d components of the Wigner 
energy from experimental data. Experimental values are consistent with the simple relation 
dT=o/W ~ 1. This is an empirical argument that the Wigner term emerges from the T=0 
part of the np pair interaction. The presence of an extra correction to the binding energy of 
odd-odd nuclei (the d-teim) may have important consequences for the independent quasi- 
particle theory of np pairing where even-even and odd-odd nuclei are treated on the same 



footing. 

To investigate the microscopic structure of the Wigner term, we performed shell-model 
calculations for even-even nuclei while attenuating, in various ways, the T=0 part of the 
effective shell-model interaction. This study confirms that the bulk part of Wigner energy 
comes from the T=0 np pairs. We also found that the contribution from deuteron-like ( J=l, 
T=0) pairs is by no means dominant; it is similar to, and near the middle of the /7/2-shell 
even weaker than, the contribution from maximally aligned pairs. The contributions from 
the T=0 pairs with intermediate spins are also non-negligible, especially in the upper part 
of the sd-shell. Finally, we found that the residual interaction between the valence proton 
and the valence neutron in an odd-odd nucleus exhibits a very weak mass dependence. 
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TABLES 

TABLE L Results of the least-squares fits for the Wigner energy strength TT(j4) = ovi//^" for 
a=l/2, 2/3, and 1. The actual minimum of the mean standard deviation a corresponds to a=0.92 
(o"~0.196). All data points with Z>10 shown in Fig. |^ were considered in the fit. 
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FIGURES 

FIG. 1. Experimental values of W (filled circles) and d [Eq. (|5|), open circles] in N=Z nuclei 
extracted from measured binding energies p^. For even-even nuclei, the values of W were obtained 
from Eq. (^), while the values for odd-odd nuclei were obtained by means of a similar indicator, 
involving binding energies of N=Z and N=Z+2 even-even cores. The triangles mark the values of 
d calculated using experimental binding energies of the lowest T=0 states in odd-odd nuclei. The 
solid line represents the average value Wav=47/A (see text). 

FIG. 2. Calculated displacement ew of the binding energy of ^^Mg and ^'^Cr from the average 
parabolic {N — Z)^ behavior along an isobaric chain. Shell-model calculations were performed 
in the Ohuj configuration space. The results of calculations for the binding energies of even-even 
nuclei along the A=A8 chain (normalized to ^^Cr) are shown in the insert. The values of ew were 
obtained using the shell-model Hamiltonian with the J = 1, 2, .., J^axj^ = matrix elements 
removed. For instance, the result for Jmax=3 corresponds to the variant of calculations in which 
all the two-body matrix elements between states \jiJ2JT=0) with J=l,2,3, were put to zero. The 
results are normalized to the full shell-model value e^y**' (</max=0). The Coulomb contribution to 
the binding energy has been disregarded. 

FIG. 3. The strength of the Wigner term, W, extracted using binding energies calculated with 
the Ohio shell model. Full shell-model calculations (filled circles) agree very well with experimental 
data (open circles). The results of shell-model calculations with the (J=l, T=0) two-body matrix 
elements removed (Jmax=l variant, triangles), with all T=0 matrix elements removed (Jmax=7 
variant, diamonds), and with (J=l,3,5,7, T=0) two-body matrix elements removed {J^^=7 vari- 
ant, squares) are also shown. In the /p-shell, only two points have been calculated due to practical 
limitations. Also a bridge between sd and fp shells cannot be covered as it requires interactions 
that include cross-shell matrix elements which are currently not established. The insert shows the 
values of W extracted from the ETFSI mass formula |2^. They are practically zero for all nuclei 
considered. 



10 



I I I I I I I I I I I I I I I I I I I I I I I I I 



> 









2.5 
2.0 
1.5 
1.0 
0.5 
0.0 






47 



(MeV) 



© 




o* 



6 



o- 



o 



• W 

o d 

^ dT=o 

•6 



A. 

d 



o 



I I I I I I I I I I I I I I I I I I I I I I I I I 



24 32 40 48 56 64 

A 



W. Satula et al., Fig. 1 






W 



1.0 

0.8 

0.6 
0.4 

0.2 

0.0 



48Cr 



25 U 




1 


1 1 

A: 


1 1 1 
=48 


1 1 


— 


V 


• 

1 1 


i 1 


1 


1 



-a- 



-o- 



/ 



-4 4 N-Z 

""•O Q - 



24Mg 



C^ 



01234567 



max 



W. Satula et al., Fig. 2 



I I I I I I I I I I I I I I 



T 



> 

9i 



^ 



2.0- 



1.5- 



1.0- 



0.5- 



.8-.«. 



§ 1.0- 



"I — I — I — I — I — r 



ETFSI 



A--A 



A' 



«:« 



0.0- 



- y:-o: 



'-^ r 



0.0 -.^^'^.•..^•^•w^-.....^ . 

I I I I I I 

40 60 80 100 

O 

•--8 

A---A 



sd fp 



o 


exp. 




• 


total 




^ ''max" 


:1 


D 


T(odd)_ 
•'max" 


:7 


~ •'max" 


:7 



I I I I I I I I I I I 'I' I I I I 



I I I I 



20 



25 



30 35 45 

A 



50 



W. Satula et al., Fig. 3 



